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ABSTRACT 

We describe techniques concerning wavelength calibration and sky subtraction to max¬ 
imise the scientific utility of data from tunable filter instruments. While we specifically 
address data from the Optical System for Imaging and low Resolution Integrated Spec¬ 
troscopy instrument (OSIRIS) on the 10.4 m Gran Telescopio Canarias telescope, our 
discussion is generalisable to data from other tunable filter instruments. A key as¬ 
pect of our methodology is a coordinate transformation to polar coordinates, which 
simplifies matters when the tunable filter data is circularly symmetric around the op¬ 
tical centre. First, we present a method for rectifying inaccuracies in the wavelength 
calibration using OH sky emission rings. Using this technique, we improve the abso¬ 
lute wavelength calibration from an accuracy of ^ 5 A to 1 A, equivalent to ^ 7% 
of our instrumental resolution, for 95% of our data. Then, we discuss a new way to 
estimate the background sky emission by median filtering in polar coordinates. This 
method suppresses contributions to the sky background from the outer envelopes of 
distant galaxies, maximising the fluxes of sources measured in the corresponding sky- 
subtracted images. We demonstrate for data tuned to a central wavelength of 7615 A 
that galaxy fluxes in the new sky-subtracted image are ^ 37% higher, versus a sky- 
subtracted image from existing methods for OSIRIS tunable filter data. 

Key words: galaxies: distances and redshifts - galaxies: clusters - astronomical 
instrumentation: interferometers 


1 INTRODUCTION 

A Fabry-Perot interferometer, or etalon, is comprised of two 
reflecting plates working in a collimated beam. For a specific 
incidence angle of incoming light, the etalon transmits light 
of wavelength A in a circular pattern of radius r around the 
optical centre. The range of wavelengths transmitted by the 
hlter is adjusted by changing the separation between the 
rehecting plates. 

Tunable hlter instruments (TFs), often built with 
Fabry-Perot interferometers, are proving to be a hexi- 
ble and cost-effective implementation of spectrophotome¬ 
try. The ability to precisely tune to an unlimited number of 
wavelengths in a specihed interval circumvents the need to 
purchase arbitrary narrow-band hlters (Bland-Hawthorn & 
Jones 1998). TFs are suitable for studies of emission and ab- 
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sorption lines in any redshift window, and they yield higher 
resolution {R ~ 500) than low-resolution grisms (Gonzalez 
et al. 2014). However, the varying wavelength across the held 
of view makes data from TF instruments challenging to deal 
with. Background sky emission can be highly variable across 
an image in which bright OH sky emission lines appear as 
prominent rings (see Section]^ for an example). Full utilisa¬ 
tion of TF data requires a precise wavelength calibration and 
robust means of subtracting the complicated sky pattern. 

In this paper, we discuss rehnements to the wavelength 
calibration and sky subtraction for TF data from the red 
mode on the Optical System for Imaging and low Reso¬ 
lution Integrated Spectroscopy instrument (OSIRIS, Cepa 
et al. 2013a,b) on the 10.4 m Gran Telescopio Canarias 
(GTC) telescope. We specihcahy consider data of emission 
line galaxies from the OSIRIS Mapping of Emission-line 
Galaxies in A901/2 (OMEGA) survey (Chies-Santos et al. 
2015) in the Space Telescope A901/2 Galaxy Evolution Sur- 
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vey (STAGES) field (Gray et al. 2009). Our discussion is 
generalisable to similar TF instruments. We summarise the 
most important properties of the data in Section Sec¬ 
tions and 1^ address the wavelength calibration and sky 
subtraction, respectively. 


2 THE OMEGA SURVEY 

Here, we briefly summarise the survey design and rele¬ 
vant data acquisition details of OMEGA. For the com¬ 
plete details, see Ghies-Santos et al. (2015). OMEGA is 
based on a 90 h ESO/GTG Large Programme allocation 
(PI: A. Aragon-Salamanca). It was designed to yield deep, 
spatially resolved emission-line images and low-resolution 
spectra covering the Ha and [Nil] lines for galaxies in the 
STAGES supercluster. 

An exposure time of 600 s was adopted to achieve an 
S/N ^ 10 in the line flux for galaxies with 17 ^ R ^ 23.5 
(Vega) continuum magnitudes. Deblending the Ha and [Nil] 
lines required a TF full width at half maximum (FWHM) 
bandwidth of 14 A and a wavelength sampling of 7 A. The 
7615-7734 A wavelength range was covered in 18 increments 
to probe the full cluster velocity range. Twenty telescope 
pointings were used to map 0.18 deg^ of the supercluster. 

Observations were taken over three observing seasons 
between 2012 and 2014. Glear sky conditions were required, 
and the moonlight was grey or dark. Different exposures for 
a given held were not necessarily observed in the same night 
or under consistent environmental conditions (e.g., temper¬ 
ature, humidity, seeing). The median seeing was ~ 079, and 
always ^ l72. Additional details concerning these observa¬ 
tions and the data reduction are provided in Ghies-Santos 
et al. (2015). 


3 WAVELENGTH SELF-CALIBRATION 

Gonzalez et al. (2014) show that the radial dependence of 
wavelength for the OSIRIS red TF is given by the expression 

A = Ao-5.04r^ + 03 (A)r'^ (1) 

where 

a 3 (A) = 6.0396 - 1.5698 x 10“®A + 1.0024 x 10“’'A^ (2) 

Ao is the effective wavelength at the optical centre (i.e., the 
wavelength to which the TF is tuned), r is measured in 
arcmin, and wavelengths are measured in A. After applying 
the above calibration to our data, we still found signihcant 
wavelength offsets between the spectra of galaxies imaged 
independently in partially overlapping helds. The magnitude 
of the offsets varied from held to held, but it was in general 
enough to ahect hux calibration and velocity measurements. 

Assuming the radial dependence of wavelength in Equa¬ 
tion (1) is correct (which we will test later in this section), 
we attempt to update the Ao term based on the positions of 
sky rings in the images. Adjusting Ao in this way essentially 
corrects for instrument tuning inaccuracies. 

The high-resolution (R 35,000 at 7000 A) spectral at¬ 
las of Osterbrock et al. (1996) shows multiple OH emission 
lines populate the spectral range of our observations. We 
therefore simulate how the sky spectrum should look given 



Wavelength (A) 


Figure 1. The dashed curve is an intermediate-resolution {R = 
2500) sky spectrum from the GTC. The solid curve results after 
convolving the dashed curve with a Gaussian kernel having a 
FWHM equal to the resolution of our data (14 A). The peaks 
of the sky lines in the convolved spectrum are measured as the 
local maximum around each peak. The colour coding and reversed 
wavelength scale make for easy comparison with Figures and 


our chosen TF bandwidth (14 A). We convolve an OSIRIS 
sky spectrum of resolution higher than our data with a 14 A 
FWHM Gaussian kernel. Figure shows the original sky 
spectrum and the result of the convolution; central wave¬ 
lengths of the sky lines in the low resolution spectrum are 
measured simply as the local maxima of the peaks. The rel¬ 
ative strengths of the night sky emission lines are known to 
vary with time, and this will affect the adopted convolved 
wavelength of the blended lines, limiting the accuracy of 
the wavelength calibrations. To evaluate the variability of 
this effect, we also convolved an independent sky spectrum 
taken with the European Southern Observatory’s Ultraviolet 
and Visual Echelle Spectrograph (UVES). All sky lines ob¬ 
tained after convolving the UVES spectrum agree to within 
±0.14 A of those from the convolved OSIRIS sky spectrum. 

We generate sky spectra for every exposure to compare 
with the sky spectrum in Figure (see Ghies-Santos et al. 
2015 for the observing strategy). For the sky background 
we simply use an intermediate-step frame from the OSIRIS 
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Offline Pipeline Software (OOPS, Ederoclite 2012) that has 
been through all processing steps (overscan subtraction, bias 
subtraction, flat fielding) except sky subtraction. (Alterna¬ 
tively, one could also directly use the sky models discussed 
in Section]^. 

The sky images are converted from Cartesian x — y coor¬ 
dinates to polar r — 0 coordinates, where r is the distance of 
a pixel to the optical centre and 0 is the angle from the im¬ 
age 2 /-axis. The conversion of Cartesian to polar coordinates 
is made by backward mapping. A grid in polar coordinates 
with the desired resolution in r — ^ is initialised, and then 
to each r — 0 pixel the intensity at the corresponding Carte¬ 
sian pixel is assigned. We have used a resolution of 1 pixel 
(0^25) in r and 1 degree in 0. Note, we adopt the optical cen¬ 
tre reported by OSIRIS handboolQof Xq = 772, Tb = 976 
(CCDl) and Aq = -35, To = 976 (CCD2). 

The example transformation in Figure shows the sky 
emission rings become vertical columns in the r — 0 plane. 
We have checked that there is no systematic change in the 
column centres (i.e., tilt) with 6. This means the sky emis¬ 
sion rings have circular symmetry and that we can accu¬ 
rately characterise their centre with a single measure after 
collapsing the r — 0 plane in 6. 

The two-dimensional images are collapsed into one¬ 
dimensional spectra by taking the median across all 0 at 
a given r. The example in Figure shows a sky spectrum 
for an image where the central wavelength is approximately 
7643 A. Comparing Figure!^ to Figurewhile noting the 
central wavelength of 7643 A, implies the visible sky lines in 
the spectrum correspond to wavelengths 7624.85, 7571.80, 
7524.75, and 7477.8 A. The central wavelength of each sky 
line in a one-dimensional spectrum is measured simply as 
the local maximum emission. 

Reliable observations in the OSIRIS red TF data are 
limited to a circular field of view of radius 4 arcmin (960 pix¬ 
els); beyond 4 arcmin, there can be contamination by other 
orders. For all sky lines at radii less than 4 arcmin from the 
optical centre, we measure the expected wavelengths of the 
sky lines using Equation (I). The average offsets between the 
actual and expected sky line wavelengths are the requisite 
adjustments to Ao for a single exposure. As a simplification, 
we calculate average adjustments to Ao as a function of tun¬ 
ing wavelength and field. The adjustments were typically 
~ 5 A, but they were as high as ~II A in some cases. 

In Figure |4(a)| we test whether the updated calibra¬ 
tion still follows the relation from Gonzalez et al. (2014, our 
Equation I). We plot Asky — Aq, where Agky is the sky line 
wavelength and Aq is the recalibrated tuning wavelength, for 
sky lines within 4 arcmin of the optical centre. The line is the 
relation from Gonzalez et al. (2014). Offsets of ~ 2 — 3 A 
from the calibration are common. Gonzalez et al. (2014) 
assumed the CS-IOO Fabry-Perot controller in OSIRIS is 
strictly linear in its gap spacing-control variable (Z) relation, 
and attributed all non-linearities to phase dispersion effects 
in the dielectric coatings of the etalons. The discrepancy we 
measure here may be an indication that the assumption of 
linearity is not strictly true. 

While the relation from Gonzalez et al. (2014) is accu¬ 
rate over a large wavelength range, systematic offsets up to 
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Figure 2. This image is a two-dimensional sky map after con¬ 
version to polar coordinates. The wavelength at the centre of the 
image is approximately 7643 A, and the actual wavelength probed 
declines from the centre outward. The vertical columns of rela¬ 
tively higher intensity correspond to sky emission rings. These 
columns show no systematic tilt and demonstrate the sky emis¬ 
sion has circular symmetry. 


Wavelength (A) 

7643 7637 7622 7594 7356 7505 



Figure 3. This one-dimensional spectrum is the result of median 
collapsing across 0 the sky map in Figure[^ The peaks correspond 
to sky lines at 7624.85, 7571.80, 7524.75, and 7477.8 A. The latter 
two lines are beyond a radius of 4 arcmin and are not used in the 
recalibration. 

3 A can occur in certain narrow wavelength ranges. One can 
proceed with this level of disagreement if it does not affect 
the accuracy of the science, e.g., Ha-based star-formation 
rates. For redshift determinations, however, a 2 A error (93 
km/s) is large. 

We propose an iterative method to refit the wavelength 
calibration and reduce the typical error down to 1 A. In our 
data set, most (89%) exposures with more than one sky line 
inside 4 arcmin of the optical centre contain the sky line 
at 7624.85 A. We focus on these exposures and perform an 
iterative procedure. 
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• Calculate offsets relative to the Gonzalez et ah (2014) 
calibration for sky line 7624.85 A as a function of tuning 
wavelength setting (Ao) and field. 

• Apply these offsets to the remaining sky lines at wave¬ 
lengths other than 7624.85 A. To these sky lines, fit the same 
functional form justified by Gonzalez et al. (2014), but with 
more free parameters for better agreement, namely 

A = Ao + a2r^ + a3(A)r^ ( 3 ) 


where 


as{\) = 03,0 - 1.5698 x + 1.0024 x 10“’'(4) 


The free parameters are a 2 and a 3 ,o, which are fixed to -5.04 
and 6.0396 by Gonzalez et al. (2014). 

• Iterate until the model converges with the sky lines 
at 7624.85 A. In subsequent iterations, the offsets for the 
7624.85 A sky line are calculated relative to the newly fit 
model in the above step. 


Figure 4(b) [ shows the result after 10 iterations of this 
procedure. The fit is noticeably better overall (although 
residuals for a few measurements slightly worsen). Most 
(85.5%) of the individual sky line measurements agree to 
within 1 A of the new calibration. This is an improvement 
over Figure |4(^ where that percentage is only 50.7% for the 
Gonzalez et al. (2014) solution. 

Using the new calibration from Figure 4(b)| we calcu¬ 
lated average corrections to Ao for each combination of wave¬ 
length setting and field in our data. Correcting the central 
wavelengths with the offsets calculated from the sky lines in 
this way yields a wavelength calibration accurate to within 
0.5 A for most (76%) of the individual image frames, and 
accurate to within 1 A for 95% of the frames. The 1 A level 
of accuracy is a success considering it is ~ 7% of the 14 A 
instrumental resolution. The calculated corrections do not 
vary in any systematic way with date of observation, am¬ 
bient temperature, or humidity. The mean offset does vary 
weakly with Aq. The average offset declines from 2 A at 
Ao = 7620 A to -2.7 A at Ao = 7734 A. 

It is important to point out that OSIRIS uses a non¬ 
standard phase correction scheme that may affect the gen- 
eralisability of this method to very different Fabry-Perot in¬ 
terferometers. Additional prerequisites for the application 
of this method are that the data be circularly symmet¬ 
ric around the optical centre and that wavelength depen¬ 
dence on detector position be radially symmetric. Our wave¬ 
length recalibration benefited from having a common sky 
line across most exposures; not having this may produce 
poorer results. This technique’s accuracy may further be 
limited by variation in sky line relative intensities, which 
can perturb the effective peak positions of the blended sky 
lines in the low-resolution OSIRIS sky spectra. For our TF 
bandpass (~ 14 A) and spectral range, the variability is 
small (±0.14 A), but it may be worse in other instances. 


4 SKY SUBTRACTION 

In this section we overview a new sky subtraction technique 
for TF data. The current sky subtraction in OOPS works by 
artificially dithering images on a 3-by-3 pixel grid (Edero- 
clite 2012). Median combining the dithered images produces 


an estimate of the sky background. Visual inspection of the 
sky maps, such as the one shown in Figure |5(a)| tuned to 
a central wavelength of 7615 A, shows that the outer en¬ 
velopes of galaxies remain in the sky background. Applying 
this sky background oversubtracts and eliminates the real 
outer envelopes of galaxies. 

We have developed an alternate approach to sky sub¬ 
traction that mitigates the problem of source oversubtrac¬ 
tion. The basic idea is to remap the flat-fielded images from 
Gartesian to polar coordinates (Section and then run a 
median filter, in polar coordinates, along circular arcs of 
grouped radii, to filter out sources and leave behind an es¬ 
timate of the sky background. 

Gircular symmetry in the background sky emission 
means that the sky changes quickly in r but is relatively 
stable across ^ at a given r (Figure]^. However, sky sub¬ 
traction is not as straightforward as converting a collapsed 
one-dimensional sky spectrum (e.g.. Figure]^ into a two- 
dimensional circularly-symmetric sky model. Because in 
some cases the sky varies as a slow function of 0 (5-10% 
amplitude) due to imperfect ffatfielding, this simplified ap¬ 
proach is less effective at removing broadened sky line emis¬ 
sion than the method we will adopt. 

The simplest application of a median filter is to filter 
over ^ at a given r with a fixed window size (as measured in 
polar image pixels). This is also not recommended. Sources 
at small distances from the optical centre are not effectively 
removed because, for a fixed window size in polar coordi¬ 
nates, the number of pixels from the original image that con¬ 
tribute to the median becomes very small at small radii. A 
significant improvement in filtering performance is achieved 
by fixing the median filter window to contain an approxi¬ 
mately constant number of pixels from the original image. 
Thus, for a given filtering window drdO^ the number of po¬ 
lar image pixels in the filter varies inversely with radius, and 
this ensures a sufficient number of pixels from the original 
image are filtered over at small radii. 

As we show below, changing the width of dO has little 
effect on the outcome of the sky subtraction. Ghanging dr, 
however, has a more significant effect. Increasing dr incor¬ 
porates more image pixels into the filter, improving robust¬ 
ness and increasing the speed of the filtering. Raising dr 
too high yields a sky map with too low spectral resolution 
that introduces artifacts in the sky-subtracted images. As 
a compromise, we have chosen to set dr to 3 pixels (0775, 
0.57 A at r = 4 arcmin), small enough to be beneficial but 
not problematic. 

Because the dependence of A on r is approximately 
quadratic, a more sophisticated approach would be to 
change dr according to quadratically varying bins such that 
each radial bin represents a constant AA decrement in wave¬ 
length. We have tested this binning pattern with our data. 
Even if we make the radial bins 1 pixel wide at large radii, 
the radial bins near the optical centre become so wide that 
often include many sources which are not effectively removed 
by the median. This leaves behind residual background light 
in the sky-subtracted images that would bias the photome¬ 
try of centrally located sources. Other, more complex, radial 
binning patterns could also be applied, but we did not find 
they yield any significant improvement over our method. 
Nevertheless, while this increased complexity is not neces- 
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Figure 4. In panel [(a)] (top), the difference in the sky line wavelength (Agky) and the recalibrated tuning wavelength (Aq) is plotted 
against radius from the optical centre. The dashed line is the wavelength calibration derived by Gonzalez et al. (2014). The bottom panel 
shows the corresponding residuals. The bands of near zero residual are from exposures with only one sky line at r < 4 arcmin. Panel [(1^ 
is similar but instead uses the result of the iterative refitting procedure discussed in Section!^ 


sary with our data, it may prove useful with different data 
sets. 


As a further refinement, we apply sigma clipping to pre¬ 
vent the extended halos of bright sources from being smeared 
into the background sky, which causes local oversubtraction 
of the sky and leads to dark halos around galaxies. We have 
found that iteratively clipping to 3cr for 5 iterations is suffi¬ 
cient to nullify this effect. 

Figure 5(b)| shows the results of applying the filtering 
to the same exposure highlighted in Figure |5(a) For this 
example, we have fixed dr to 3 pixels and the total filter size 
drdO to 170 pixel-radians, corresponding to a filter size of 
20 degrees at r = 2 arcmin. The sky background is much 


smoother than what is provided by the current algorithm in 
OOPS (Figure [5(a)[ left-hand panel), and the sky-subtracted 
sources in the right-hand panel of Figure |5(b)| are clearly 
brighter than in Figure |5(a)| Binning in r, however, yields 
sky models with band-like artifacts of width dr (see Figure]^ 
row 2, column 2). The change in brightness between bands 
is typically ~l-2% of the sky background, smaller than the 
background shot noise (3-4%). 

Figure gives a magnified view of the sky subtraction 
in an image subregion containing bright extended objects as 
well as small faint sources. The sky subtracted-images in the 
third column emphasise that galaxy light sacrificed by the 
standard OOPS algorithm is retained with median filtering. 
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(b) 

Figure 5. For a representative exposure, panel [(a)] shows the sky background measured by OOPS (left), and the corre spon ding sky- 
subtracted image (right). The wavelength at the optical centre is approximately 7615 A. For the same exposure, panel shows an 
example of the sky map generated by the filtering technique outlined in Section (left) and the corresponding sky-subtracted image 
(right). Median filtering was performed with dr set to 3 pixels (0775), a total filter size drdO of 170 pixel-radians, and iterative Scr 
clipping. The colour bars show the colour stretch and are similar in both panels. 
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Figure 6. A magnified view of F igures [5 (a) | and |5 (b) | showing a subregion containing both bright extended sources and compact faint 
sources. The top row shows the original flat-fielded image still containing the sky, the sky background from OOPS, and the final OOPS 
sky subtraction. The second row shows the sky background inferred from Altering and the resulting sky-subtracted image. The colour 
bar shows the colour stretch and applies to all panels. 


To quantify how much light was gained with the re¬ 
vised sky subtraction procedure, we measured the fluxes of 
sources in the sky-subtracted images using the methodology 
of Chies-Santos et al. (2015). Fluxes of galaxies in the right- 
hand panel of Figure [5(b) are higher by a median of ~ 37% 
versus the OOPS image in Figure |5(a)[ This boost in flux 
is not strongly sensitive to either dr or the total Alter size 
drdO. We repeated this test for Alters with dr values of 1-5 
pixels (at Axed drdO of 170 pixel-radians), as well as for drdO 
values of 10-30 degrees at r = 2 arcmin (with dr = 3 pixels). 
The resulting median flux ratios deviate by ~ 1 — 3% from 
the ~ 37% value obtained with the adopted hducial Alter 
parameters {dr = 3 pixels, drdO = 170 pixel-radians). 

This Altering technique is implemented in Python using 
standard libraries, including Astropy (Astropy Collabora¬ 
tion et al. 2013). We are assisting with the incorporation of 
this technique into OOPS so that ah users of OSIRIS will 
have access (Ederoclite, private communication). 


5 SUMMARY 

In this letter, we have used data from the red TF mode on 
OSIRIS to demonstrate new techniques for wavelength cali¬ 
bration and sky subtraction. Central to our methodology is 
the use of polar coordinates, which simplihes matters when 
the TF data is circularly symmetric. In Section we out¬ 
lined a technique for wavelength recalibration using OH sky 
emission rings. This approach increases the accuracy of the 
absolute wavelength calibration from ~ 5 A to 1 A, or ~ 7% 
of the instrumental resolution. In Section we presented 
a new method to estimate the sky background by median 
Altering in polar coordinates. The merit of this approach is 
that of light from the extended halos of emission-line galax¬ 
ies does not contaminate the background sky maps. Sources 
in the associated sky-subtracted images will likely be signih- 


cantly brighter than how they would appear in images based 
on the current sky-subtraction algorithm in the OSIRIS 
reduction software pipeline, OOPS. Future OSIRIS/OOPS 
users will beneht from this sky subtraction method. 
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